Vortex formation in two dimensions: When symmetry breaks, how big are the pieces? 
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We investigate the dynamics of second order phase transitions in two dimensions, breaking a gauged 
U(l) symmetry. Using numerical simulations, we show that the density of topological defects formed 

when the dynamics is overdamped at the instant 



-1/2 



scales with the quench timescale tq as n -~ iq 

when the freezeout of thermal fluctuations takes place, and n 
This is predicted by the scenario proposed by one of us El. 
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in the underdamped case. 



PACS numbers: 05.70.Fh, 11.15.Ex, 67.40.Vs, 74.60-w 



Phase transitions occur at all energy scales, from Bose- 
Einstein condensation near absolute zero to the sub- 
Planck temperatures relevant for symmetry breaking in 
a cosmological setting. While the dynamics of some 
types of first-order transitions (i.e., the process of nu- 
cleation) has been extensively studied and is well under- 
stood for some time, analogous understanding of second- 
order transitions is emerging only now. In this case there 
is no supercooling and the final state of the system is 
asymptotically approached through phase ordering. Un- 
till recently, research has focused largely on the asymp- 
totic scalings of this post-transition ordering rather 
than on the dynamics of the transition itself. 

The change in focus is relatively recent. Kibble, in 
a seminal paper |||, pointed out that topological de- 
fects may have significant cosmological consequences - 
i.e. they may act as seeds for structure or as constraints 
on models. While their distribution at the time they ap- 
pear is largely forgotten by the present, some of its fea- 
tures (such as the ab initio existence of 'infinite' string) 
are essential for the scenario. Moroever, the initial den- 
sity of topological defects may be directly relevant for 
generation of baryons j| . 

The initial configuration of the order parameter estab- 
lished in the course of the transition is therefore impor- 
tant. Furthermore, it is accessible in cases where topo- 
logical defects are formed, as they bear witness to the 
dynamics of the order parameter in the immediate vicin- 
ity of the critical point. Experiments based on this idea 
allow one to probe the critical dynamics of symmetry 
breaking, and have been carried out in liquid crystals || 
and in superfluids ■ This has already led to new in- 
sights into the dynamics of the transition between phases 
A and B of 3 He 0. 

We study the consequences of second order phase 
transformation of the order parameter tp, a complex 
scalar field, with Landau- Ginzburg dynamics in two spa- 
tial dimensions and the associated gauge field A a . This 
is the Abelian Higgs model coupled to a heat bath and 
with a dissipative term. When cast into first order form, 
the equations of motion are the following (tpx — C01) "02 )j 
i runs over {x, y}, and D a ip = d a + ieA a ip); 



kx = dti>x, n 4 = d t Ai 

d t TT X = VVx - e 2 A 2 ip x - 2ee XY A l d l ^ Y 
dV 

Wx + v 

d t n t = V 2 A, - eexvipxdiipY - e 2 A^\ 2 
V(v» = -l-em& 2 + ]iP 4 . 



(1) 

(2) 
(3) 
(4) 



The system is subjected to the white noise D(x,y,t); 
(d{x,y,t)d{x',y',t')) = 2r,6S(x - x')S(y - y')S(t - t>). 
Here 9 is the heat bath temperature and r\ sets the rate at 
which the field is damped, in accord with the fluctuation- 
dissipation theorem. 

We induce the symmetry breaking by changing the 
sign of the dimensionless parameter e in the effective po- 
tential over the quench timescale tq, so that e = t/rg 
(N < !)• The phase transition occurs when it be- 
comes energetically and entropically favourable for the 
order parameter to assume (in equilibrium) a finite ex- 
pectation value. In our case this happens in the region 
< £ « 1. The shift of critical temperature from e = 
occurs as a consequence of the coupling to the gauge field 



-e 2 (A 2 }) and as a result of finite temperature 



9, which we take to be 0.01. 

The equations of motion are evolved on a square lattice 
of 512 2 grid points, using the standard staggered leapfrog 
method with periodic boundary conditions. We impose 
the gauge A t — 0, implicit in Eqns. (2) and (3), and aver- 
age each case over 20 realisations. Defects were resolved 
by several grid spacings at low temperatures. 

We begin the simulations well above the phase transi- 
tion, allowing the system to come to equilibrium under 
the influence of the noise and relaxation at constant e. 
We monitor the behaviour of the order parameter and 
the gauge field throughout the subsequent evolution. The 
focus of attention, however, is the number of topological 
defects, which can be identified as zeroes of i\) in the bro- 
ken symmetry phase. Initially, in the symmetric phase, 
such zeroes are plentiful and short lived (see Fig. 1). 
While they cannot be identified with defects, their den- 
sity and arrangement gives an idea of the nature of the 
thermal fluctuations. 
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(a) (b) (c) 

FIG. 1. Representations of the field and defects (filled squares or open triangles) on 128 2 grid during and after a 
quench of timescale tq = 32, with 77 = 1.0 and in the presence of a gauge field (e = 0.5). Figure (a) shows the field 
at the critical point (e = 0), (b) at e = 0.6 and (c) at e = 1, time t = 2tq after the critical point. 



Below the transition their overall density decreases, 
although there still exist regions where the field has a 
near-zero expectation value and, hence, plenty of unsta- 
ble zeros. Eventually, only a few isolated defects (as well 
as a few pairs which are about to annihilate) remain. 
We count only zeroes that have no companions within 
Co = mo" 1 . 

The local and global (e = 0) gauge cases are quali- 
tatively indistinguishable until this late stage. However, 
local defects do not interact significantly over more than 
a few correlation lengths and annihilate relatively slowly, 
even when the friction coefficient is small. By contrast, 
global defects interact via a potential logarithmic in dis- 
tance, and disappear much more rapidly. Estimates of 
initial defect densities become more difficult in this case. 
Below, we shall focus primarily on the local case, leaving 
the detailed global-local comparisons for the future. 

The theory of defect formation combines the realisa- 
tion, due to Kibble ||, that the domains of the order 
parameter ip which break symmetry incoherently must 
contain of the order one 'fragment' of a defect, with the 
estimate jTJ of the relevant size based on the comparison 
of the relaxation timescale of ip with the effective rate 
of change of the mass parameter e. In the immediate 
vicinity of the critical temperature, the dynamics of ip 
are subject to critical slowing down. Thus, the timescale 
r over which the order parameter can react is given by 



TQ 

|e|V2 



(r = l/m ) (5) 



in the overdamped and underdamped cases respectively, 
where correspondingly the first or second time deriva- 
tive in Eqn. ^ dominates. The overdamped scenario 
is presumably more relevant for condensed matter ap- 
plications, while in the cosmological settings ip may be 
underdamped. 

The characteristic timescale of variations of e is e/e = t. 



The system is able to readjust to the new equilibrium as 
long as the relaxation time is smaller than t. Hence, 
outside the time interval [— t, t] defined by the equation 
r(e(t)) = t, the evolution off/' is approximately adiabatic, 
and physical quantities follow their equilibrium expecta- 
tion values. Thus the times 



and 



-2/3 1/3 
T Q 1 



H = ±1 



± 
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(6) 



(7) 



mark the borders between the (approximately) adiabatic 
and impulse (or, perhaps, "drift") stages of evolution of tp 
in the overdamped and underdamped cases, respectively. 

In particular, the correlation length £ of if) above 
the transition will cease to increase with the Landau- 
Ginzburg scaling (£ = Co/kl 1 ^ 2 ) above the transition 
once the adiabatic-impulse boundary at —t is reached. 
Dynamics will be suspended (except for the drift and 
noise) in the interval [— t, t] and will resume at +t below 
the transition, as the mass term drives the symmetry 
breaking. 

We expect, then, that the characteristic length scale 
over which ip is ordered already in the course of the 
transition will be the correlation length at freeze-out, 
£ = £o/Ve 0- This results in 



tq 
^t 2 



1/4 



^ - Co 



1/3 



(8) 



in the two cases. The initial density of vortex lines in 
two dimensions should then scale as 



1 



1 



(9) 
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(10) 



in the over- and under-damped regimes respectively. 
Above, / is the proportionality factor which is expected 
to be of the order of a few ||, and which may be esti- 
mated analytically in some cases pT| . 

The relative importance of the ip and %j> terms in Eqn. 
(||) also depends on e. What matters for the formation of 
defects is - in view of the arguments above - which of the 
two terms dominates at t. Thus, n 2 > e is the condition 
for overdamped evolution, which leads to the inequality 
rfrQ > uiq. Therefore, one can enter the overdamped 
regime by performing a sufficiently slow quench, as well 
as increasing the damping parameter n. 

We verify these scalings in Fig. 2, where the numbers 
of defects obtained in both the overdamped and under- 
damped cases are plotted function of the quench 
time tq. In the overdamped case the annihilation of the 
well-separated topological defects is slow. Consequently, 
it is relatively easy to count them at some fixed time. 
This yields Fig. 2b with n ~ t%, 7 = -0.44 ± 0.1 for 
T Q > 8, in good agreement with the theoretical predic- 
tion of 7 = — 1/2 [l] . A similar conclusion can be reached 
for the underdamped case, Fig. 2a, except that annihila- 
tion is now quite rapid, and the initial number of defects 
is harder to define. This is especially true for the fastest 
quenches which produce most defects. However, when 
the two left-most points most affected by annihilation 
are ignored, straightforward counting of vortices yields 
7 = —0.6 ± 0.07. This is in agreement with the theory - 
which in this case predicts 7 = — 2/3 [JL 12| - and with the 
indications from kink formation in one dimension | fL3| . 

The reason the slope may be expected to become more 
shallow for small tq is easy to understand. A very fast 
quench becomes indistinguishable from an instantaneous 
one, which does not allow for the adiabatic regime we 
have noted above. Instantaneous quenches start the evo- 
lution in the broken symmetry phase, but with an ini- 
tial field configuration which will contain many zeros 
per healing-length size volume - too many to regard 
them as well-defined defects. Unless the defects are well- 
separated at their conception (which in effect requires 
6«1), annihilation will decide their initial density. We 
see this already in Fig. 2b, where e = 0.5 for the left-most 
point. The effect is even more pronounced in the under- 
damped case, which allows for more rapid annihilation 
(Fig. 2a). 

The annihilation rate can be expected to be of the 
simple form n = —xn 2 , where x is a function of n, as the 
annihilation is proportional to the rate at which defects 
encounter one another. This leads to 



X t. 



(11) 



We can then infer the value of the initial defect density 
no, which is also shown in Fig. 2a. This procedure (fol- 
lowed, for instance, in the superfluid work, JjJ) yields a 
somewhat different slope no ~ r -o.79±o.o4^ s ^ ee p er than 
equal-time data, and somewhat steeper than the theoret- 
ical prediction. It also seems to succesfully correct for the 
annihilation (although the left-most point still seems to 
be affected). Also, using the data in Fig. 2, we estimate 



/ ; ~ 12 and / , 



(equal-time measurements), / r ~ 4 



(from fitting to Eqn. [Tl] ) . 

One of the issues in the formation of topological de- 
fects is the relative significance of the thermal fluctua- 
tions, which continue to re-arrange the configuration of 
the order parameter down to the Ginzburg regime (i.e. 
below the phase transition temperature). If, as was once 
thought, fluctuations and Ginzburg temperature were to 
determine initial defect density, then taking the system 
above the critical point would be expected to erase the 
pre-existing configuration of defects and create a new 
one. We have performed a numerical experiment de- 
signed to test the importance of fluctuations. The initial 
configuration with defects present is the one shown above 
in Fig. lc. We reheat this by varying e from the broken 
symmetry value e = 1 to several values in the vicinity of 
the critical point e ~ 0, using the same tq was used as 
in the original quench of Fig. 1 . 

The results (Fig. 4) demonstrate that even when the 
system is taken above the critical point, the initial config- 
uration of defects is eventually recovered, as long at the 
re-heating does not take it further than e into the sym- 
metric phase. This confirms the theory put forward in 
Q, and leads one to conclude that thermal fluctuations 
cannot significantly re- arrange configurations of the order 
parameter on scales larger than ~ £, unless the "impulse 
strip" |e| < |e| is traversed (or unless the time spent in 
that regime is set by a timescale other than the original 
quench timescale tq). 

We have used high-resolution numerical simulations to 
explore phase transitions in two dimensions and have 
found the scaling with quench timescale and damping 
agree with the predictions of the Kibble-Zurek scenario. 
The importance of the freeze-out time t as the defining 
moment for defect formation has been illustrated. 

The authors would like to thank Pablo Laguna for the 
use of his code and Mike Warren for assistance with the 
Loki parallel computing facility at Los Alamos. 



Fig. 3 shows (in the underdamped case, when the anni- 
hilation is appreciable) the fit between (111]) and the data. 



[1] W.H. Zurek, Nature 317, 505 (1985); Acta Phys. Pol. 

B24, 1301 (1993) 
[2] e.g. A. Bray, Adv. Phys. 43, 357, (1994) 
[3] T.W.B. Kibble, J. Phys.: Math. Gen. A9, 1387 (1976); 



3 



Phys. Rep 67, 183 (1980) 

W. Perkins et al, Nucl. Phys. B353, 237 (1991) 
P. Laguna & W.H. Zurek, Phys Rev. Lett. 78, 
(1997); see also cond-mat/9705141 



2519 



I. Chuang et al, Science 251, 1336 (1991); M. J. Bowick 
et al, Science 263, 943 (1994) 
P.C. Hendry et al, Nature 368, 315 (1994) 
V.M.H. Ruutu et al, Nature 382, 332 (1996); 
V.M.H, Ruutu et al, Phys. Rev. Lett., submitted; 
cond-mat/970603S| (1997) 



C. Baiierle et al, Nature 382, 334 (1996) 
Yu. M. Bunkov, O. D. Timofeevskaya, 



cond-mat/9706004 



A.J. Gill, G. Karra, R. Rivers, in 
Lythe, Phys. Rev. E53, 4217 (1996) 
W.H. Zurek, Phys. Rep. 276 177 (1996) 
P. Laguna and W. H. Zurek, |hep-ph/9711411 



preparation; G.D. 



4 




100 1000 10 100 



FIG. 2 (left) The variation of defect 
count N with tq in both (a) under- and 
(b) over-damped regimes. The fitted 
slopes are (a) -0.60 ± 0.07 (\ 2 = 0.034, 
dropping the 3 left-most points) and (b) 
-0.44 ± 0.10 (x 2 = 0.44, dropping 2 
points) Predicted values were —2/3 and 
— 1/2. The points fitted by the dashed 
line are inferred from the fits in fig. 3. 
The slope is -0.79 ± 0.04 ( x 2 = 0.27). 
FIG. 3 (right) Defect annihilation in the 
underdamped regime. The inverse of the 
defect count is plotted against time for 
various tq. Also shown are least-squares 
fits using Eqn. 11, 
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FIG. 4. Memory and reheating. The post-transition configuration shown in 
figure 1 (the point S in the sketch to the left) is reheated to various temperatures 
close to the critical point (figs, a-c) and cooled again to e = 1 (A-C respectively). 
The reheat values of e are (a) -0.1, (b) -0.2 and (c) -0.25. Here |e| ~ 0.17. 
Notice that the memory of the configuration is largely preserved even when the 
critical temperature is exceeded during reheating (figs. A and B). Memory is 
erased only when the temperature crosses the |e| freezeout zone associated with 
the formation of the original defect configuration. 
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